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Abstract 



Multiscale stochastic volatility models have been developed as an efficient way to capture the principle 
effects on derivative pricing and portfolio optimization of randomly varying volatility. The recent book 
Fouque, Papanicolaou, Sircar and S0lna (2011, CUP) analyzes models in which the volatility of the 
underlying is driven by two diffusions - one fast mean-reverting and one slow-varying, and provides 
a first order approximation for European option prices and for the implied volatility surface, which is 
calibrated to market data. Here, we present the full second order asymptotics, which are considerably 
more complicated due to a terminal layer near the option expiration time. We find that, to second order, 
the implied volatility approximation depends quadratically on log-moneyness, capturing the convexity 
of the implied volatility curve seen in data. We introduce a new probabilistic approach to the terminal 
layer analysis needed for the derivation of the second order singular perturbation term, and calibrate to 
S&P 500 options data. 

1 Introduction 

Stochastic volatility models relax the constant volatility assumption of the Black-Scholes model for option 
pricing by allowing volatility to fluctuate randomly. As a result, they are able to capture some of the well- 
known features of the implied volatility surface, such as the volatility smile and skew. While some single- 
factor diffusion stochastic volatility models such as Heston's [T3], enjoy wide success due to the existence of 
semi-analytic pricing formula for European options, it is known that such models are not adequate to match 
implied volatility levels across all strikes and maturities; see, for instance, Numerous empirical studies 
have identified at least a fast time scale in stock price volatility on the order of days, as well as a slow scale on 
the order of months, for example [31 [51 [13 HZ] • This has motivated the development of multiscale stochastic 
volatility models, in which instantaneous volatility levels are controlled by multiple driving factors running 
on different time scales. 

A class of multiscale stochastic volatility models is analyzed in [7] , where an approximation for European 
options and their induced implied volatilities is derived, which can capture the overall level of implied 
volatility, its skew across strike prices and its term-structure over a wide range of maturities. However, the 
analysis there is limited to a first order approximation, which cannot pick up the slight convexity of the 
observed equity implied volatility surface. In this paper we extend the results of [7] to second order. This 
extension is non-trivial, as it requires a careful terminal layer analysis, which we approach probabilistically. 
For some related multiscale perturbation techniques in European option pricing, we refer for instance to [3] 
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and [1] (spectral methods), [TB] (matched asymptotic expansions), pQ, [T3] and [TU] (Malliavian calculus), 
and [21] (inner-outer expansions). 

Our second order results allow us to capture the slight convexity of the implied volatility skew. Ad- 
ditionally, we are able to maintain analytic tractability which is important for calibration to data, as we 
demonstrate. Of course, numerous asymptotic regimes have been analyzed in recent years for the option 
pricing problem in incomplete markets: see [5], [T2] and [TH] for some references. Here our focus is not just 
on deriving and proving convergence of the approximation in the appropriate limits, but in disentangling the 
calibration procedure that results from it. Compared to the first order theory, this is much more involved 
as there are many more group parameters and basis functions that have to be accommodated to implied 
volatility data. Despite the increase in complexity, we show this can be implemented successfully. 

The rest of this paper proceeds as follows. In Section [31 we describe the class of multiscale stochastic 
volatility models that we will work with. Using a formal singular and regular perturbation analysis, we 
derive a pricing approximation which is valid for any European-style option. We establish the accuracy of 
our pricing approximation in Theorem 12.41 In Section [31 we present an explicit formula for the implied 
volatility surface induced by our option pricing approximation. Additionally, we show how a parameter 
reduction, crucial for calibration purpose, can be achieved with no loss of accuracy. In Section 13.21 we 
outline a procedure for calibrating the class of multiscale stochastic volatility models to the empirically 
observed implied volatility surface of liquid calls and puts. We carry out this calibration procedure on call 
and put data taken from the S&P500 index. Section @] concludes. 



2 Second Order Option Pricing Asymptotics 

We consider the class of multiscale stochastic volatility models studied in [8]. Let X denote the price of 
a non-dividend-paying asset whose dynamics under the historical probability measure P is defined by the 
following system of stochastic differential equations (SDEs): 



dX t = pX t dt + f(Yt,Z t ) X t dW x , 
dY t = -a(Y t ) dt+-^= f3(Y t ) dW v tl 

dZ t = S c(Z t ) dt + Vd g(Z t ) dW t z . t 



Here, (W x , W v , W z ) are P-Brownian motions with correlation structure 

d{W x ,W y ) t = p xy dt, d(W x ,W z ) t =p xz dt, d{W\W z ) t = Pyz dt, 



where (p xy , p xz , p yz ) satisfy \p xy \,\p xz \,\p yz \ < 1 and l + 2p xy p xz p yz - p 2 xy - p 2 xz - p 2 yz > 0, which guarantees 
that the correlation matrix of the Brownian motions is positive-semidefinite. The asset X has geometric 
growth rate p and stochastic volatility f(Y t , Z t ) which is driven by two factors, Y and Z. Under the physical 
measure, the infinitesimal generators of Y and Z are scaled by factors of 1/e and 8 respectively. Thus, 
e > and 1/6 > represent the intrinsic time-scales of these processes. We will work in the regime where 
e << 1 and 6 << 1 so that Y and Z represent fast- and slow- varying factors of volatility respectively. Most 
importantly, we assume the fast factor is mean-reverting, that is that Y is an ergodic process with a unique 
invariant distribution n under P, which is independent of e. 

Under the risk-neutral pricing measure P (chosen by the market) the dynamics are described by 



dX t =rX t dt + f{Y t ,Z t ) X t dW x , 
1 



dY t 
dZ t 



-a{Y t ) - -^=A(Yt, Z t ) p(Y t ) ) dt + -L f3(Y t ) dW t 



(Sc(Z t )->/8T(Yt,Zt)g(Ztj) dt + Vs g(Z t ) dW z , 



(2.1) 



2 



where (W x , W v , W z ) are P-Brownian motions with the same correlation structure as between their P- 
counterparts, and r is the risk-free rate of interest. The functions A(j/, z) and T(y, z) represent market prices 
of volatility risk, which we have assumed such as to preserve the Markov structure of (X, Y, Z) and of (Y, Z) 
by itself. 

Consider a European option with expiration date T and payoff }i(Xt)- The no-arbitrage price of this 
option at time t < T can be expressed as a discounted expectation of the option payoff 

P £ >*(t, X t , Y t , Z t ) = 1 [e- r(T -^h(X T )\x t , Y t , Z, 

Here, E denotes an expectation taken under the pricing measure P. Note that we have used the Markov 
property of (X, Y, Z) to express the price of the option as a function P £ ' S (t, x, y, z) of the current time t and 
the state variables (X t ,Y t , Z t ). Using the Feynman-Kac formula, one finds that P e ' S satisfies the following 
partial differential equation (PDE) and terminal condition: 

L**P s ' s = 0, P e > s (T,x,y,z) = h(x), (2.2) 

where, introducing the notation 

T> k = x k d*... x , fc=l,2,..., (2.3) 

the operator L £,s is given by 

L e ' s = (^L + -t=£i + L^j + y/S (Af^a + +5M 2 , 

with 



£ = -p 2 {y)d 2 vy + a{y)d y: 



1 

2' 

<Ci = p xy f3(y)f(y, z)Vid y - /3(y)A(y, z)d y , 

L2 = d t + y 2 (y,z)'D 2 +rV 1 -r, 

M 3 = Pyz p{y)g{z)d 2 yz , 

Mi = p xz g(z)f(y, z)D x d z - g(z)T(y, z)d z , 

M 2 = l -g 2 {z)dl z +c{z)d z . 

For general (/, a, /?, A, c, g, T) no explicit solution to (|2.2p exists and we seek an asymptotic approximation 
for the option price. The fast factor asymptotics is a singular perturbation problem, while the slow factor 
expansion is a regular perturbation. Thus, the small-£ and small-5 regime gives rise to a combined singular- 
regular perturbation about the 0(1) operator L 2 - We expand P e,s in powers of y/e and y/S as follows 

P £ > s (t, x,y,z) = ^2J2 yft'yft'Pijb ^> y, z). 

j>Q i>0 

This is a formal series expansion, and our approach is to find P{ t j for i + j < 2 and establish an accuracy 
result for the truncated series. We point out that we are working within an infinite-dimensional family 
of models since the functions (/, A, T) are unspecified: the 18 group parameters that are found in Section 
12.71 and calibrated in Section 13.21 contain specific moments of these functions identified by the asymptotic 
analysis. 

Because we are performing a dual expansion in half-integer powers of e and 5, we must decide which of 
these parameters we will expand in first. We choose to perform a regular perturbation expansion with respect 
to 6 first. Then, within each of the equations that result from the regular perturbation analysis, we will 
perform a singular perturbation expansion with respect to e. As the combined regular-singular perturbation 
expansion is quite lengthy, to aid the reader, we provide a summary of the key results in equation (|2.41[) . 
located at the end of this section. 
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2.1 Regular Perturbation Analysis 

The regular perturbation expansion proceeds by expanding L 6,5 and P e,s in powers of VS as follows 



L £ ' s =L e + VSM £ + SM 2 , P e =5^V£ J 'Pj, (2.4) 

j>0 



where from 



£ £ = -£ + ^£i+£ 2 , M £ = ^M 3 + M 1; Pf = VVeVij. (2.5) 

e Je Je 4 



i>0 



Inserting (|2 .4[) into (|2.2p and collecting terms of like-powers of Vs, we find that the lowest order equations 
of the regular perturbation expansion arc 



0(1) 
0(6) 



= £ £ P £ , (2.6) 
= £ e Pf + M £ P £ , (2.7) 
= £ £ P 2 £ + M £ Pf + M 2 P £ . (2.8) 



Within each of these three equations, we now perform a singular perturbation analysis with respect to e. 
2.2 Singular Perturbation Analysis of the 0(1) Equation ( 12. 6j) 

Beginning with the 0(1) equation, we insert expansions (|2.5|) into (|2.6|) and collect terms of like-powers of 
J~e. The resulting 0(l/e) and 0(1/ Je) equations are: 

0(l/e): = £ P ofi , 

0(1/ J^): = £ o Pl,o+£iPo,o. 

We note that all terms in £o and <&i take derivatives with respect to y. Thus, if we choose Pq^o and Pi,o 
to be independent of y, the above equations will automatically be satisfied. Hence, we seek solutions of the 
form 

Pa.o = Po,o(t,x,z), Pifi = Pi,o(t,x,z), 

i.e., no y-dependence. Continuing the asymptotic analysis, the 0(1), 0(y/e) and 0(e) equations are: 



0(1) 

0(V£) 

0(e) 



= L P 2 ,a +U&Tfi+ ^Pxo, (2.9) 
= £0^3,0+^1^2,0 + £2^1,0, (2.10) 
= L O P 4 , O + £iP 3 ,o + £2^2,0, (2.11) 



where we have used the fact that L\P\fi = 0. 

Equations (|2.9j) . ()2.10|) and (|2.11j) are Poisson equations of the form 

= L o P + X- (2.12) 

By the Fredholm alternative, equation (|2.12|) admits a solution P only if \ is in the orthogonal complement 
of the null space of the adjoint operator Lq. Since the unique invariant distribution II satisfies £qTI = 0, 
this leads to the solvability or centering condition: 

(x) ■= I x(y)II(dy) = 0. (2.13) 
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Applying the centering conditions to equations 
P\fl do not depend on y, we find 

0(1) 
0(e) 

where the operator (£ 2 ), given by 



(|2.10[) and (|2.11[) . and using the fact that Po,o and 



= (£ 2 ) Po,o, 

= ^iP 2> o) + (£2) Pi,o, 

= (£1^3,0) + (£2^2,0} , 



(£2) = d t 



\d 2 (z)'D 2 + r'D 1 -r, 



(2.14) 
(2.15) 
(2.16) 



(2.17) 



with 



a\z) :=(f(.,z)) = f 2 (y,z)n(dy) 



We observe that (£2) is the Black-Scholes pricing operator with effective averaged volatility tr(z) = y/ a 2 (z) 
in which the level z of the slow factor appears as a parameter. 

Expanding the terminal condition in (|2.2[) leads to the terminal conditions 



0(1) 



fo.oC^Xjz) = /i(x), 
Pi,o(T,x,z) =0. 



(2.18) 
(2.19) 



We denote the solution to (|2.14l) with terminal condition (|2.18p by 

Po,o(t, x, z) = P B s(t, X] a(z)), 

the Black-Scholes price of the option with volatility a(z), maturity T, and payoff function h. 

In order to make use of equation (|2.15[> , we need an expression for (£>iP%fi). Using (|2.14[> . we re-write 
(|2~9)) as follows 

£ P 2 ,o = -£ 2 Po,o = - (£2 - (£2)) Po,o = -\ (f - (f 2 )) D 2 Po,o. 



Introducing a solution 4>(y, z) to the Poisson equation 

£ o = / 2 -(/ 2 ), 

we deduce the following expression for P 2 $: 



(2.20) 



P2,o{t,x,y,z) = --^(y,z)T> 2 P ,o(t,x,z) + ^2,o(*, x, z), 



where F 2i a(t,x, z) is independent of y. 



(2.21) 



Remark. The form of (|2.21[) shows that the natural terminal condition P 2 q(T, x, y, z) = is not enforceable 
because the singular perturbation with respect to the fast factor creates a terminal layer near t — T . However, 
as we will demonstrate in Section 12 .61 the ergodic theorem enables us to impose the averaged terminal 
condition 



0(e): 



(P 2 , (T,X,-,z)) = 0. 



(2.22) 



and to obtain the desired accuracy of our pricing approximation. In fact, we will see that this is the only 
appropriate choice for proof of convergence. 
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Remark. The solution of the Poisson equation (|2.20p is defined up to a constant in y. We choose this 
constant by imposing the condition 

(<K;z))=0, (2.23) 
and we will show in Section 12.61 that this choice does not affect the accuracy of our pricing approximation. 
Inserting (|2~2Tj) into ([2~T5]) yields the following PDE for P ifi 

(£3) Pi,o = - (£iP 2 ,o) = - / (pxyPfVrfy -pkd^J (~<f>V 2 Pofi + F 2 , ^ = -VP ,o, (2.24) 

where the z-dependent operator V is given by 

V(z) = V 3 («)DiO) 2 + V 2 («)D 2) 

and we introduce the notation 
1 



V 2 {z) = - (/3(-)A(-, Z)) , V 3 (Z) = -~ Pxy (P(-)f(; Z)d y 4,(; z)) . 

The solution P10 of the PDE (|2.24l) with terminal condition (|2.19l) will be given in Proposition 12. 31 

To determine P 2i q, given by (I2.21[) . we need a PDE and terminal condition for the unknown function 
F 20 . The terminal condition F 2i0 (T, x, z) = is imposed by averaging (12.211) . and using (|2.22p and (|2.23[) . 

Proposition 2.1. The function F 2i o(t, x, z) satisfies the following PDE and terminal condition 

(£ 2 )F 2 .o = -AP Qfi -VPi,o, F 2fi (T,x,z) = 0, (2.25) 

where the z- dependent operator A is given by 

A(z) = A 2 (z)V 2 1 V 2 + A 1 {z)T) 1 V 2 + A (z)D 2 + A(z)D 2 2 , 

A 2 {z)= l -pl y {P{-)f{;z)d v M-,z)), 

M*) = ~\p*v «J3(-)H;*)8vM',*)) + (#•)/(•. *)dvM',z))) > 

A {z) = \{P{-)K{;z)d y i> 2 {;z)), 

A{z) :=-\{ (4>(; z)f(-,z)) - (<f>(;z)) (f(;z)) ) , 

and ipi(y, z) and tpi(y, z) satisfy 

£0 fa = Pfdy^ - (0fd y cf>) , L V-2 = Pkd y 4> - (pAByfi , (2.26) 

Proof. The proof is in appendix [XJ □ 

The solution F 2 ,o of the PDE with terminal condition (|2.25j) will be given in Proposition 12.31 This is as 
far as we will take the asymptotic analysis of the 0(1) equation 



2.3 Singular Perturbation Analysis of the 0(y5) Equation ( 12.71) 

Proceeding as in Section l2~2l we insert expansions (|2.5p into (|2.7p and collect term of like-powers of y/e. The 
resulting 0(VS/e) and Q(y/5/y/e) equations are: 

0(VS/e): = £ o Pq,i, 
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where we have used M^Po^ = since M 3 contains d y and Po,o is independent of y. Recalling that all terms 
in Lq and L\ also contain d y , we seek solutions Po,i an d P1.1 of the form 

Pq,i = Po,i(t,x,z), Pi,i = Pi t i(t,x,z). 

Continuing the asymptotic analysis, the Q(V~5) and 0(\/Sy/s) equations are: 

0(VS): = £ P 2 ,i +i^Pra+ £2^0,1 +^Pifi + M 1 P 0fi , (2.27) 

O(VS^): = L Q P 3tl +L 1 P 2A +L 2 P 1A +M 3 P 2 . + M 1 P lfi - (2.28) 

Equations (|2.27[) and (|2 . 28[) are Poisson equations of the form (|2.12p . Applying the centering condition (|2 . 13[) 
to (j2~?7l) and (f2T2"5l) yields 

O(VS): 0=(£ 2 )Po,i + (M 1 )P 0| o, (2.29) 

0(V*Vi) : = (£iP 2 ,i) + (L 2 ) Pi,i + (M 3 P 2i0 ) + (Mi) P lj0 . (2.30) 

We also impose the following terminal conditions 

0(VS): P ,i(T,x,z) = 0, (2.31) 

0(VSy/e) : P 1A (T,x,z) = 0. (2.32) 

The PDE (|2 .29[) and terminal condition f|2.31[) can be used to find an expression for Po,i, which will be given 
in Proposition 12.31 

The operator (Mi) appearing in (|2.29[) is given by 

(Mi) - Pxz g (/) 1>id z - g (V) d z = ^{V 1 (z)'D 1 d z + V (z)d z ), 

a 

where a' = d z a and we introduce the notation 

Vi(z) = \p xz a'{z)g{z) (/(•, z)) , V (z) = -\o> '(z)g(z) (T(;z)) . 

In order to make use of equation (|2.30|) . we need expressions for (£iP2,i) and (M. 3 P2 t o)- 
Proposition 2.2. We have the following expressions: 

(£iP 2 ,i) = (VsV^i + V 2 V 2 ) P 0) i + j t {C 2 D\ + C1D1 + Co) d z P ,o, (2.33) 

(M 3 P 2 , ) = ^CT> 2 d z P , , (2.34) 
a' 

where 

C 2 (z) = -p xy p xz d % '(z)g(z) (f3(-)f(-,z)d y ip 3 (-,z)) , 

Ct{z) = p xz a'(z)g(z) (P{-)kd y ^) + p xy g (/?(■)/(•, z)O y M; «)) , 

C (z) = -a'(z)g(z) (/3(-)A(-, z)8 y M;z)) , 

C{z) = -\ Pyz a'{z)g{z) (p(-)d v <t>(; z)) , 
and ip 3 (y, z) and ipiiVi z ) satisfy 

£0^3 = / - (/) , £0^4 = r - (r) . (2.35) 
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Proof. The proof is in Appendix [B] □ 
Inserting (f^33l) and (p04]) into ([2~30]> , we find 

(£2) P1.1 = -VPo.i - ^ea z P ,o - (Mi) Pl i0 , (2.36) 
er 

where the z-dependent operator C is given by 

G(z) = C 2 (z)T>\ + Ci(«)Di + C (z) + C(z)V 2 . 

The solution Pi j of the PDE (|2.36l) with terminal condition (|2.32p will be given in Proposition 12. 31 . This is 
as far as we will take the asymptotic analysis of equation (|2.7p 



2.4 Singular Perturbation Analysis of 0(5) Equation ( 12. 8j) 

We now move on to the 0(8) equation (|2.8I) . Proceeding as in Sections 12.21 and 12.31 we insert expansions 
(|2.5[) into (|2.8I) and collect term of like-powers of y/e. The resulting 0(8 /e) and 0(<5/v / e) equations are: 

0(S/e): = L P a, 

0(8/ Ve) ■ = £ Pi.2 + ^iPo,2 + jtfaPoT, 

where we have used 3Yf.3Po,i = since M3 contains d y and Po,i is independent of y. Recalling that all terms 
in £q and L\ also contain d y , we seek solutions Po i2 and Pi 2 of the form 

Po,2 = Po,2(t,X,z), Pl,2 = Pl,2(t,X,z). 

Continuing the asymptotic analysis, the 0(8) equation is: 

0(8): = £oP2,2+^Pr^+£ 2 Po.2+3tf 3 Prr+MiPo,i+M 2 Po, . (2.37) 

Equation (|2.3T|) is a Poisson equation of the form (|2.12l) whose centering condition (|2 . 13[) is 

0(8): 0-(£ 2 )P ,2 + (M 1 )Po,i+M 2 P ,o. (2.38) 

We also impose the following terminal condition 

0(8): P 0<2 (T,x,z)=0, (2.39) 

The solution P 02 of the PDE (|2.38[) with terminal condition (|2.39[) will be given in Proposition 12.31 This is 
as far as we will take the combined singular- regular perturbation analysis. 

2.5 Review of Asymptotic Analysis and Pricing Formulas 

In the previous sections we showed (formally) that the price of a European option can be approximated by 

(2.40) 

where 



pe,S 




Po,o + Ve Pi,o + VS Pq,i 


+ eP 2 ,o + 


8P 2 , Q + Vl8P 1A 




0(1): 


(£ 2 ) Po,o 


= 0, 




P ,o(T,x,z) = 


h(x)," 


Q(VS): 


(£2) Pi,o 


= -VPo.o, 




P x ,o(T,x,z) = 


0, 


0(V8) : 


(£ 2 )Po,i 


= - (Mi) P ,o, 




Po,i(T,x,z) = 


0, 


0(e): 


P2,0 


= -i</)2) 2 P ,o + P 2 ,o, 






> 




(£2) P2,0 


= -yiP ,o-VPi,o, 




F 2 . (T,x,z) = 


0, 


0(8): 


(£2) Po,2 


= - (Mi) P ,i - M 2 P , , 




P . 2 (T,x,z) = 


0, 


0(Ve~8) : 


(^2> Pl,l 


= -VP ,i-^e9 2 P ,o- 

it' 


(Mi) Pi i0 


, P 1A (T,x,z) = 





(2.41) 
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and the ^-dependent operators in (|2.41[) are given by 

1 



0(1): 




0(Vi): 


V 


0(VS) : 


(Mi) 


0(e): 


yi 


0(6): 


M 2 


0(V7S) : 


e 



2 

i 7 



(2.42) 



In the following, we provide explicit expressions for {Pi.j, i + j < 2}. 

Proposition 2.3. Introducing r = T — t, we have the following expressions for the {Pi,j}: 
Po,o = P B s{cr(z)), P h0 = TVP BS (a(z)), P ,i =T7fid tr P B s{v(z)) 



"2,0 

Po,2 
Pl.l 



- cj)D 2 PBs(o-(z)) + F 2 ,o, w/iere P 2 ,o 

2r 2 - 2 
3 



+-r^ )P BS (a(z)), 



^KiKi 9 CT + + + T -B 2 (% a + ±=d}j + T -B X dj\ Pbs(°(z)), 

(r 2 VN 1 d r , + Leda + ^XiVj Pbs(«(z)), 



where we have introduced the z- dependent operators 

Ni = Vi Di + V , K = V{ Di + Vq, 

and parameters 



V; = d z V 3 , j = 0,1, 2, 3 
Proof. The proof is in Appendix [C] 



s 2 = ~<?V) 2 , 



V = Vg'DiDa + F 2 'D 2 , 



Bi = ~<?V'W. 



□ 



2.6 Accuracy of the Approximation 

To establish the accuracy of our pricing approximation P e ' S defined in (12.401) . we make the following as- 
sumptions: 

1. The system of SDEs (|2.1I) has a unique strong solution (X, Y, Z) for fixed e,5 < 1. 

2. The market prices of volatility risk are bounded: ||A||oo < 00 and ||r||oo < 00. 



3. Let be a diffusion process whose infinitesimal generator is Lq (so that, in distribution, Y t — Y$ 

under P). We assume that is ergodic and has a unique invariant distribution II with density n, 
and that £0 has a positive spectral gap. We note that two of the processes that are most commonly 
used to model volatility — the Cox-Ingersoll-Ross (CIR) and Ornstein-Uhlenbeck (OU) processes — 
satisfy these assumptions. 

4. The process Y^' admits moments of any order uniformly bounded in t: 

kl 



sup E 

t 



Y, 



(i) 



< C(k). 
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5. Let Z^> be a diffusion process whose infinitesimal generator is M2 (so that, in distribution, Z t = Z 
under P). We assume that Z^ admits moments of any order uniformly bounded in t < T: 



(i) 
st 



supE 

t<T 



(1) k 



<C(T,k). 



6. We assume that the function f(y, z) is smooth in z, that the solutions <^(y, z) and {ip%{y, z), i < 4} to 
equations (I2.20|) . (|2.26j) and (|2.35|) . are at most polynomially growing in y and z, and that 



a\z) := / f'(y,z)U(dy) < 00. 

7. We assume that the payoff function h(x) is smooth and bounded with bounded derivatives. The proof 
of accuracy provided here (in Appendix [D| uses this smoothness assumption. It can be generalized to 
nonsmooth payoff functions (such as in the case of call options considered in Section [3]) following the 
regularization argument given in [6] for the first order approximation with a fast factor. The details of 
the generalization to the second order approximation with fast and slow factors are quite lengthy and 
we omit them here. 

Theorem 2.4. For fixed t < T , x, y, and z, the model price P £ - 5 solution of fl2H) and our price approxi- 
mation P e ' S defined by (|2.40[) satisfy 

\P E ' S (t, x, y, z) - P E ' S (t, x, y,z)\ — 0(e 3/2 + eV5 + 8y/£), 
Proof. The proof is given in Appendix [D] □ 

Remark (Terminal Layer Analysis). The main difficulty in extending the accuracy of our pricing approxi- 
mation from first order to second order is the treatment of the terminal condition for the second order term 
P2,o arising from the singular expansion due to the fast factor Y. In |16) . the solution i-2,0 is derived by a 
formal matched asymptotic expansion with a terminal layer of size e. Here, in Appendix [D] we provide a 
probabilistic proof based on the ergodic property of the fast factor Y, which justifies the choice of terminal 
condition made in (|2.22l) . 

2.7 Group Parameters and Exotic Option Pricing 

We now summarize the parameters needed in the pricing approximation formulas derived in the previous 
section. We begin by separating the y-dependent part in P e ' S given by (I2.40p . by writing 

P e >\t, x,y,z) = ~e cf>(y, z)V 2 P , {t 7 x, z) + Q^ s (t, x, z), 

where 

Q £ ' 5 (t, x, z) := P ,o + Vei\,o + ^ Pas + V76 P M + e P 2 , + S P , 2 . (2.43) 

Using (|2.40p . (|2.41[) and the linearity of the operator (£2), we find that Q £ ' S satisfies the following PDE and 
terminal condition 

(L 2 ) Q e ' s = S e ' s , Q £ < 6 (T, x, z) = h(x), 

where the source term S e,d is given by 



s £ < s = - Vivp ,o - Vs (Mr) p ,o - VTs (vp ,i + ^ed z p ,o + (Mi) Pi, 

- e {A P .o + V P lj0 ) - S ( (Mi) P ,! + M 2 P ,o) 

= - (7iV)P ,o - (VS (Mi))P ,o - (\/iV)(V£P .i) - (V^C) ±jd z P . - (VS (Mi))(Vi-Pi,o) 

a 

- {eA)P , - (ViV)(ViPi,o) - (VS {Mi))(VS P 0A ) - {6M 2 )P , Q . 
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To extract which group parameters are needed for the price expansion, we absorb a half-integer power of s 
and/or S into the corresponding group parameters and define: 

V?:=y/iV t , V*:=V5Vi, A\:=eA u B\ := S B u Cf :=V^Ci. (2.44) 

Similarly, we absorb the appropriate e or <5 pre-multiplier into the terms of the expansion (|2.43l) by defining 
Pf and Pi, through 

yfeP lfi {t,x,z) = P[ >0 (t,x;a(z),Vi(z),Vi(z)), VdP ,i(t, x, z) = P^x; a(z), V s (z), V?(z)). 

Substituting from (|2.42p the expressions for M 2 , V, A, (Mi) and C, and changing the d z derivatives in (Mi) 
and M2 acting on Po,o into d a derivatives acting on Pbs(ct(z)), we finally have 

S s ' s = - (ViV^ + V 2 £ V 2 ) P BS - 2 (V?D a + V S ) d a P BS 

- {y^ x V 2 + V £, D 2 ) P* tl - (c £ 2 ' Sr D\ + CfV x + C £ ' S + C E ' 5 D 2 ) d a P BS 



- 2 (yfDi + V S ) U + ^fd Vf + %dv^j Pto 

- (AiD?D 2 + Al-D^z + A £ Q V 2 + A £ vfj P BS - (Vf Di© 2 + V 2 £ V 2 ) Pf : 




Here our notation is V i ' £ (z) — d z Vf(z), and similarly V( 5 . Since P^o is linear in V3 and V 2 and P0.1 is linear 
in Vi and Vb, neither dv^Pi.o, dv?Pi,o, d v sPo,i nor dysP^^i contain any of the V^s (that is, they are order 
one quantities). 

As such, the group parameters that appear in the source term S £ ' S and therefore, in the price approxi- 
mation (|2.40p are 

V 3 £ , V 2 £ , Vi, V S , Cf, cf, Cf, C £ > s , A £ 2 , AIA £ ,A £ ,B S 2 ,B(, YLYLYlL^ YL. (2 .45) 

g o~ a a 

These 18 parameters, which move with the slow volatility factor Z t , as well as <jf{y,z) := e<j){y,z) needed 
in (|2.40j) . can be obtained by calibrating the class of multiscale stochastic volatility models to the implied 
volatility surface of (liquid) European options, as described the Section 13.21 Note from (I2.44[) that the Vf 
are order y/e, the Vf order \/~5 and that they appeared in the first order asymptotic theory in .7 . The 
new parameters (A £ , Bf, C £ ' ) come from the order e, order 5 and order \/eS terms in the the second order 
expansion respectively. 



2.7.1 Parameter Reduction 

The group parameters in (|2.45|) depend on the current level z of the slow volatility factor and, in the case of 
(jf , on the fast factor too. In order to calibrate completely from the implied volatility surface and not use 
historical returns data to estimate a(z), we replace it by a quantity cr*(z) which absorbs the term V 2 (z). In 
so doing, there is now one less parameter, and we show that the accuracy of the second order approximation 
is unchanged. The accuracy of our approximation is given in Theorem 12.41 and the following Proposition 
refers to it. 

Proposition 2.5 (Parameter Reduction). Without loss of accuracy, in the pricing equations and formulas 
obtained in this section, we can replace a(z) by o~*(z) defined by 

a*(z) := y/a(zf + 2Vi(z), (2.46) 

replace p^o by the Black-Scholes price at volatility cr*(z), and remove the V 2 -dependent terms. 

Proof. The proof is in appendix U3 □ 
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2.7.2 Pricing Exotic Options 



It is important to remember that the formulas given for the Pij's in Theorem 12.31 are valid for European 
options only. This is because Theorem 12 .31 relies on the Vega-Gamma relationship (|C.1|) . which does not hold 
in general for exotic options. Nevertheless, the PDEs derived in (|2.4ip are valid for many path-dependent 
options with appropriate boundary conditions (e.g. barrier options). The approximate price of a path- 
dependent option can by found by numerically solving an inhomogeneous PDE (with additional boundary 
conditions) , the main point being that the parameters involved are the same as those appearing in European 
option prices used for calibration as described in Section 13.21 A full development for particular contracts is 
outside the scope of this paper. 



3 Asymptotics for Implied Volatilities and Calibration 



In what follows, we give the second order expansion of Black-Scholes implied volatility, which is obtained by 
inverting the Black-Scholes formula for European call options with respect to the volatility parameter. 

3.1 Implied Volatility Expansion 

In practice, it is common to find the parameters of a pricing model by calibrating the model to the implied 
volatility surface of liquid European calls and puts, rather than by calibrating to call and put prices directly. 
To this end, we seek an implied volatility expansion 



j>0 i>0 



1,3 



such that 



P £ ' 5 = Pbs (I 5 ' 5 ) 



Performing a Taylor expansion of Pbs(P' 5 ) about I and rearranging terms yields 

Po,o + y/e Pi,o + P ,i + VeSPi,i + e P 2 ,o + 6 P , 2 + ■ ■ ■ 

= Pbs{Iq,o + \fehfl + VSI ,i + VeS h,i + eI 2 .o +8Io,2 H ) 

= Pbs(Io,o) + fve/i,o + V8I ,i +■ ■ ■) d a P BS (I ,o) 



I (V^Ii,o + V5I OA + --) 2 d? (r P BS (I , 



= Pbs(Io,o) + y/ihjad a P B s{hja) +^ I ,xd a P BS {h) 
+ VFS (h, Q I 0A dl a P BS (I Qfi ) + h&PBsilofi)] 



^O.l^o-er-P-BsUo.o) + ^0,2^0- P B s(^0,f 



Equating terms in p. II) of like powers of \fe and y/5, and using Pq^q = P B s{o~) we find 
0(1) : J„,o - a, 0(e) : h,o 



P 2 ,o 1 T 2 dl&Po,o 

o 1,0" 



0(Ve): h,o = 



Pifi 
daPo,o ' 
Po 



0(8) : Jo,2 



9crPo,0 2 9 a Pofi 

Pay _ 1 j2 dgo-Po,o 
d<jPo,o 2 0,1 d a Poo 



daPa.o 



daPoj 



daPofi 



(3.1) 



(3.2) 
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For a European call or put option with strike price K and time to maturity r it is convenient to express the 
-fjj's as functions of forward log-moneyness 



d := log (K/xe rT ) (forward log-moneyness). 

Setting the payoff function h{x) = (x — K) + for a call option and using the expressions given for {Pij} in 
Theorem 12. 31 the Iij's in (I3.2[) become 

0(1) : 7 0>0 = a, 



0(Vi): h,o = V 2 ~ + V 3 (1= 



a \2a ■ 
t d 



0{VS): Io,i = Vot + Vi(^ 



a 



2 r a \ 2 a 6 J \ r a° 2 a 6 

3d 



Vo 



A 2 



f 3d 2 


3 


\ t 2 <j 7 


+ 2ra 5 2 


' d 2 


1 - M 


t 2 a 5 


t a 3 4 a J 


( ^ ~ 


1 d 


\t 2 a 5 


T (T 3 T<7 3 



" : : 1 1?)-* 



1 



0(5) : / ,2 = V? — + VbVi -— T + — + V 2 ' 



6 a 7 \3er 3 6cr/ \ 3 er 5 6cr 3 6cr 3 6 a 



M3.3) 



^oV 2r M t/ /t 3 2dr\ Vq / r 2 2d 



+ V °t - T+3^ + ^#V3 ' 3, 2 



! r 



0(V7j) : /!,! = y ^2 (-js) + V V 3 ( — rj 



V/ /t 2 2d 2 2r 2dr\ / d 2 

1 I I 1 I I o I 

3a 2 3, 
3d r 



U + 3^-3^+3^ +B2 U^ + 6^-T2- |+Sl 



r 



CT 4 2ct 

3d r \ „„ / 6d 2 3 3d 



a 4 2c? 2 y ■ v ™ 6 1 ct 4 fj 4 



a' \<jJ a' \a A 2 a J a' \a A 2 a 

( — - — + — + —\ +r ( T - + A L,J_ 

V Ira 5 a 3 a 3 4a 2 \8 2ra 4 2a 2 2a 2 



T 



d \ „ /t\ „ / t d 2 1 



4 ' 2 cr 2 / V 2 / V 8 2r o- 4 2 a 2 

Observe that this second order expansion produces an implied volatility curve which is quadratic in 
log-moneyness d and therefore accounts for the slight turn in the skew that is most prominent in shorter 
maturity options data, as we will see in Figure [TJ The first order approximation derived in [7] is linear 
in d and therefore only accounted for the skew effect. Note also that the parameter reduction outlined in 
Proposition 12.51 can be applied to this implied volatility expansion as well (a replaced by a* and Wterms 
removed) , and this will be used in the calibration in the next section. We also remark that the formal second 
order expansion for the case of a single slow volatility factor had previously been considered in [9] , |18) and 
[20] . for instance. 
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3.2 Calibration 



In this section we discuss how the parameters (|2.45l) . which are needed to price exotic options as discussed 
at the end of Section 12.71 can be obtained by calibrating the multiscale class of models to liquid European 
option data. We define 



Using 



where 



I £ ' S ■■= Io,o 



\fehfl + VS Io,i + Vs~5 + e/2,0 + S Iq.2- 
and the parameter reduction described in Proposition 12.51 we have 



Je,8 = 



0(l/r) 
0(1) 

0(r) 

0(r 2 ) 

0(d/r) 
0(d) 

0{dr) 

0(d 2 /T 2 ) 

0(d 2 /r) 
0(d 2 ) 



1 \ d d 2 

— k + I + rni + t 2 n) H — (p + Tq + T 2 s)^ ^ (u + t v + t w) , 

t I t x ' t z ' 



?V2 



3(^0 
2(cr*) 5 



1 



(a*) 3 

•ye ,5 



(ct*) 3 2ct* 

Qe,S 



{a*) 4 2(a*) 2 


2((7*) 2 






4(7* 4ct* (a*) 3 a*' 




H H + 

(7* 2(7* 


h 2a* 


#1 . °0 | °1 






2 2 4 


' 8 8 ' 6(ct*) 3 





V s Vi 



2V 1 S V{ S 



2(er*) 2 6a 
(V S ) 2 , W , iV*) 2 



3(a*) 2 a*' 



+ 



V 5 V' £ 
2a*a*' 



4a*a* 



B\a* 2V 5 V " 5 



6a* 



6a* 



q = 



Vq Vj_ 

3(7*' 

m? , 

2((T*) 5 

(a*) 4 



6cr* 

3a*' 
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6(7*' ' 



3(7*' 



1 



(a*f 



+ 



1 , 



(a*) 3! 



e,<5 



C 



e,5 



((7*) 4 2((7*) 2 2((7*) 5 



v s vi e VfV£ e Vf 



(a*) 3 a*' (a*) 3 a*' (a*)" 2 



W£V{ 5(V{ 

3((7*) 3 



(5\2 



W 



[a*y 

' (a*f 
7(Vf) 2 



+ 



6(a*Y 
(a*) 5 



2V^ S V 1 5 
3(a*) 2 a*' 
A e 



2V$V{ S 



2VfV{ 



3((7*) 2 (7*' 3((7*) 2 CT* 



r *\5 : 



+ 



s,6 



2((7*) 4 2((7*) 4 ((7*) 5 (7* 



2V{V{ S 



(3.4) 



3((7*) 5 3((7*) 3 3((7*) 4 (7*'' 

In total, we have ten "basis functions" with which to fit the empirically observed implied volatility surface: 

1 d d 2 d 2 

-, 1) t, t 2 , -,d, dr 7 — , — , d 2 

T T T T 
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It will be helpful to define 

:= {k, I, m, n,p, q, s, u, v, w}, 

<^ := {a\Vi ,V^V^Cl'\C['\Cl'\C^\Al,A\,Al,A^BlBl%,%,\,^}. 

g a a 

We let I[t, d) be the implied volatility of a European call option with time to maturity r and forward log- 
moneyness d as observed from option prices on the market. Wc let I £ ' S (t, d; O) be the implied volatility of 
a European call as calculated using p.4[) . The calibration procedure consists of the following steps: 

1. Find 0* such that 

m J n E E ( / ( T "^) - ?' < Wi;©)) 2 = EE - ^''(Ti.di; e*)) 2 , 



where the double sum runs over all maturities n and strikes Kj (corresponding to forward log- 
moneyness dj) for which a call or put is liquidly traded. This is the least-squares fit of formula 
(|3.4[) resulting in estimated k,l,m, • ■ ■ ,w. 



2. Next the ten constraints of equation (|3.5j) are used to find the minimal L 2 set of parameters <£>*. That 
is, we find <&* such that 

min ||$|| 2 = ||$*|| 2 , J = {$ : equation holds with = 0*}. 

We emphasize that our calibration procedure encompasses all maturities, that is we do not fit maturity- by- 
maturity. 

3.3 Data 

We perform the described calibration procedure on European call and put options on the S&P500 index on 
two separate dates, one pre-crisis on October 19, 2006, and one post-crisis on March 18, 2010. In Figure [1] 
we plot the implied volatility fit from October 19, 2006. The parameters obtained from the above calibration 
procedure are 

cr* = 0.2051, Vi = -0.0034, Vf = 0.0023, F 6 = -0.0064, C e / = -0.0073, C{' 6 = -0.0171, 



Co- = 0.0183, C M = 0.0047, A% = -0.0002, A\ = 0.0038, A e = -0.0183, A e = 0.0011, 

y/S y/S y/S 

B s 2 = 0.0080, B\ = 0.0183, -2- = 0.0146, ^- = -0.3104, -2-= 0.9856, £ = -0.0181. 

a' a' a' 



In Figure [2] we plot the implied volatility fit from March 18, 2010. The parameters obtained from the above 
calibration procedure are 

cr* = 0.2269, Vi = -0.0062, Vf = -0.0026, = 0.0208, C e 2 s = -0.0031, C{' S = -.00034, 

Cq S = -0.0035, C e - 5 = 0.0033, A\ = 0.0034, A\ = 0.0034, A% = -0.0004, A e = -0.0012, 

y/e y/S y/S 

B 5 2 = 0.0012, B{ = -0.0035, -4r = -0.1590, 4r = 0.0914, -2- = -0.0729, (f> £ = -0.0443. 

a a' a' 

Notice that, in both cases, the obtained parameters other than cr* are small, as expected in the regime of 
validity of our expansion (i.e., small e and small 5). 
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4 Concluding Remarks 



We have derived a second order asymptotic approximation for European options under multiscale stochastic 
volatility models with fast and slow factors. Proof of convergence requires a terminal layer analysis that 
is developed probabilistically, in contrast to the techniques of matched asymptotic expansions that are 
more common in fluid mechanics. The price approximation is translated to an implied volatility surface 
approximation which is quadratic in log-moneyness and highly nontrivial in the term structure direction. 

In applying asymptotic approximations, two questions are crucial: how is it calibrated to data, and how 
can the recovered (partial) information be used. Here we have shown that the complicated second order 
formulas fit the data well across strikes and maturities (Figures Q] and [2]). Moreover, the extracted parameters 
are small when they should be small in the regime of the asymptotic analysis (Section l3.3[) . These parameters 
can then be used to value exotic options such as barrier options, that are characterized by the same PDE 
but with different boundary conditions, to the same order of accuracy, as discussed in Section [^771 

It remains to investigate stability of the calibrated parameters over time, that is, to show empirically all 
but one of the parameters vary like the slow volatility factor, and the remaining one more rapidly. This 
detailed data work is in preparation. 



A Proof of Proposition 12.11 

The goal is to compute (£2) -^2,0- Starting from the expression (|2.21[) for P2,o, applying the operator £2 and 
averaging, we obtain: 



(£ 2 P 2 ,o) = ^£2 I --^D 2 Pq,o + F 2fi j ) = ~2 W 1 ^) ^2^0,0 + (£2) F 2fi 
= ~\ ({</>) (dt -r) + \ (<pf) D 2 J ©2^0,0 + (£2) ^2,0 

= -~ [{4>) (dt-r+^f 2 )^ +i((0/ 2 )-(0)(/ 2 ))B 2 ) » 2 Po,o + (£2) ^2,0 
= -^J^fThpZ,- i( (0/ 2 ) - {4>) </ 2 ) )b 2 P , + (£2) F 2fi 

= AD 2 2 P , + (£ 2 )P 2 ,o, (A.l) 

where A(z) is given in Proposition 12. II 

On the other hand, by (|2~T6)) . (£2^2,0) = - (£1^3,0), and therefore we compute P 3j0 . From (|2~T0|) . (|2~T5l) . 
(|2.20p . (|2.21|) . and the definitions of £1 and £ 2 , we express £o-P3,o as follows 

£o-P3,o = — (£i-P2,o + <£ 2 -Pl,o) 

= - (£iP 2 ,o - (£l^2,0» - (£2 - (£2)) Pl,0 

= -£i \ - ^©2^0,0+^2,0 



£i r-i^D 2 p 0!0 +F 2 ,o)\- Q (/ 2 - </ 2 » D 2J Pi,oJ 
= - (-\p* y (pfd y 4> - (pfdyt))®!^ + ^(jSAdyj - (pAdy^Vz) P , - (~£ o 



Using the functions tpi and tp 2 defined in (|2.26[) . the solution to the previous equation is given by 

^3,0 = \p*y ^lB 1 2) 2 P ,0 - ^22)2^0,0 - ^0B 2 Pl,O + *3,0, 
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where Fafi(t, x, z) is independent of y. Now, we can compute (£iP3,o) : 

= \ply (PfWl) ^lWo,0 - \p X y (Phdylfr) £l© 2 P ,0 - \p X y (Pfdylh) »lD 2 Po,0 

+ ~ (fiAdyih) - \px V (Pfd y 4>) DiDa^o + 5 (/3A9^) D a Pi,o 

= (a 2 B2© 2 + A 1 D 1 B 2 + A D 2 ) P ,o + (v&DiDa + V 2 D 2 ) P 1)0) (A.2) 

where A 2 (z), Ai(z) and ^4o(z) are given in Proposition l2.il Inserting (jA.ip and (|A.2I) into (|2.16p yields the 
following PDE for P 2 , 

(£2) P 2 ,o = — -A Po,o — 'VPi.Oj 

which agrees with (|2.25l) . Lastly, the terminal condition F 2 ^(T, x, z) = is imposed by averaging (|2.21l) . 
and using ([2T2"2"j) and (|2~25)l : 

(P 2 ,o(T, x, •, z)) = -ij^rD 2 P ,o(T, 1, z) + F 2 , (T, x, z) = P 2 , (T, x, z) = 0. 

B Proof of Proposition 12.21 

In order to compute (£iP 2j i) we first compute P 2 ,i- Using (|2.27[) . (|2.29p and the definition (|2.35p of ip3 and 
■04, we find 

£oP 2 ,i = — £ 2 Po,i — ^VtiPo.o 

= - (<C 2 - (£ 2 )) Po,i - (Mr - (Mi)) P ,o 

= -\ (f - (f)) 2) 2 Po,i - A«S (/ - (/}) Did z P 0fi + g(T- (T))d z P ,o 
= -^£o0I ) 2 Po,i - Px Z g^oip3^>id z Po,o + gL ip 4 d z P ,Q. 
Thus, P 2 .i is given by 

P 2 ,i = -~0D 2 P„,i ~ Pxz9^id z P 0fi + g^d z P Qfi + F 2 ,i(t,x, z), (B.l) 
where P 2i i(i,x, z) does not depend on y. Next, using expression (jB.ip for P 2) i we find 

(£iPm) = Up xy PfDi-Ph)d y ui^D a fl),i 



d z Po,o 



-\p*y (Pfd v 4>) 2)iD 2 P ,i + ^ (/3A9^) 2) 2 Po,i 

(/3fd v ip3) D\d z Pos) + Px Z g (f3Ad y ip 3 ) Di9^P 0j o 
+ Arj,5 {Pfdyi>i) I>i<9 2 Po,o - g (f}Adyip4) d z P ,o 

(y 3 2)iD 2 + y 2 D 2 ) Po.i + ^7 (C 2 2)? + + Co) d z P , , 
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which agrees with (|2.33[) . Next, using expression (|2.21[) for P2,o we find 

(M 3 P 2 , > - ((p V zP(y)g(z)d 2 yz ) (~(f>D 2 P , +F 2 , 

a' 

which agrees with (|2.34[) . 

C Proof of Proposition 12.31 

As previously noted, the operator (£2), given by (|2 . 1 T[) . is the Black-Scholes pricing operator Lbs, where the 
Black-Scholes volatility a has been replaced by an average level of volatility a, and so that Po,o = Pbs{^{z)). 
In order to derive expressions for the higher order {Pij, 1 < i + j < 2}, we need the following two lemmas. 
Recall also the notation r = T — t. 

Lemma C.l. For European options, 

d a P BS {v) = TaV 2 P BS (a). (C.l) 

Proof. This is the classical relation between Vega and Gamma for plain vanilla European options. It is easily 
obtained by differentiating the Black-Scholes PDE with respect to a to obtain a Black-Scholes PDE with 
source for the Vega which in turn can be solved explicitly in terms of the Gamma. □ 

Using Lemma fC. II and the fact that the logarithmic derivative operators T>k in (|2.3[) commute (2\.2) m = 
f m Djfc), which implies that (£2) and any commute ((£2) Dfe = T>k (£2)), one can show: 

Lemma C.2. For European options, and positive integers k and n, 

T n+1 

(£2} — -P({V k })PBs{.<r(z)) = -T n P({<D k })P BS (<t(z)), 

71+1 
T n+1 

(£2) — d a Pbs{o{z)) = -r" P({2) fc }) d a P BS (v(z)), 
n + I 

T n+1 / 1 \ 

(£2) —-P({T> k }) + — - —d a P BS (v(z)) = -r" P({V k }) dl a Pbs{o{z)), 
n + 6 \ a (n + 2) J 

where P({T>k}) is some polynomial ofDi,T>2, • • ■ ,Dfc- 

Proof. The proof is a straightforward computation. In showing the second and third relations, the d a partial 
derivatives acting on Pbs are first converted into D2 using Lemma IC.ll which now commute with any 
operators and (£2)- The final step uses that (£2) Pbs(&(z)) = 0. □ 

Using Lemmas I C . 1 1 and I C . 2 1 a direct computation shows that the {Pi.j} of Theorem l2. 31 satisfy the PDEs 
of (12.41|) and the appropriate terminal conditions. 

D Proof of the Accuracy Theorem 12.41 

In what follows, we will make use several times of the fact that Po,o and its derivatives D k Po,o are bounded 
in x (because of our smoothness and boundedness assumption on the payoff function h and its derivatives). 
We will also use the fact that Y and Z have moments of all orders uniformly bounded in e and 5 (thanks 
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to the assumptions made on and in section [2~B| . The proof of this can be found following Lemma 
4.1 of 0. 

We begin the proof by recalling our price approximation P £ ' S from (|2.40[) : 

pe,6 ^ pe,6 = p^ + ^ p^ Q + y/fi p QA + VTS P 1A + E P 2 , + <5P ,2, 

where {Pij,i +j < 2} are given in Theorem 12.31 Since the singular perturbation argument involves terms 
with higher order in e, we introduce 

pe,6 = pe,5 + £ 3/2p + £ 2p + p + £ 3/2^ p 



and we observe that 



Next, we define the residual 



P 



e,5 r>€,S 



0{e 



3/2 



R e 



pe.8 pe.S 



and therefore, the proof consists of showing that R E ' 5 = 0(e 3 / 2 + eyS + 5\fe). 

From the equations (|2T2]) , ([2^9]) . (|2T0|) . (|2~TT|) . (I^27|) . (|2T28|) . and (j2T37)) satisfied by P e ' S and < 
4, j < 2, i+j < 4}, we deduce that the residual i? ■ satisfies the following PDE: 

L £ ' 5 R e < s + 0{e 3/2 + eVS + 5y/e~) = 0. (D.l) 
From the terminal conditions for {Pij,i -\- j < 2}, we deduce the terminal condition for the residual: 



R S,S (T, x, y, z) = -eP 2 , (T, x, y, z) + 0(e 3/2 + eVtf + <Vi), 



(D.2) 



where it is important to note that the non- vanishing terminal condition P 2 fi(T, x, y, z) plays a particular role 
at the e order. The probabilistic representation of R E,b , solution to the Cauchy problem (|D.lj) - (|D.2j) . is 

R £ ' s (t, x, y, z) = -e E tiX>y>z [e~ r{T -^P 2fi (T, X T , Y T , Z T )] + 0(e 3/2 + eVS + Sy/e), 

where Et, x ,y,z denotes expectation under the (e, 5)-dependent dynamics ()2.1j) starting at time t < T from 
[x, y, z). From the expression i-2,0 = — \<i>'^2Po,o + i*2,o in (|2.21j) and the terminal condition F 2i q(T, x, z) = 
in ()2.25|) . we deduce 



R e ' 5 (t,x,y lZ ) = -E t , x , y , z 



e -r{T-t)^Y T , Zt)T> 2 Po,o(T, X t , Zt)] + 0(e 3/2 + eVS + 5y/e) 



As 8 — > 0, the process Z converges to the constant process z, and as e — > 0, the process X converges 
in distribution to a geometric Brownian motion X with volatility a(z) and independent of the Y process. 
Therefore, up to an error term of order 0(£ 3 / 2 + e\/~8 + Sy/e), we deduce that 



R e ' 6 (t, x, y, z) = e CEt, v [4>(Yt,z)] + 0(e 3 / 2 + eVS + Sy/i), 



(D.3) 



where the constant 



C := E t , x 



does not depend on (e,S). 

If A = 0, that is, if the process Y had infinitesimal generator Lq under the pricing measure P, then due 
to the ergodic theorem, Et, y [4>(Yr, z)] would converge to zero exponentially fast as e — > since ((f)) = 0. 



21 



This is exactly where we see that our choice of terminal condition for P 2j o was necessary because if ((f)) =/= 0, 
then the residual would be of order e. However, under the pricing measure P, the process Y does not have 
generator £ due to the presence of the possibly nonzero market price of volatility risk A(y, z) and therefore, 
we need to analyze the behavior as e — > of M tjV [4>(Yt, z)] for A 7^ 0. 

Without loss of generality we can assume t — 0. By rescaling time it is enough to consider E [4>(Y T / e , z)] 
under the autonomous dynamics 

dY t = (a{Y t )-^~ep{Y t )HY u z))dt + !3{Y t )dWf, Yq = y. (D.4) 

The process Y (ID.4[) admits an invariant distribution Tl e with density 

J £ / fy a(u)-^(u)A(u,z) , 

where J £ is a normalization factor. 

From the ergodic property of Y and the assumption of a positive spectral gap, for e small enough, one 
can find positive constants C\ and A such that 

|E[0(y T/£ ,z)]-(0) e |<c ie - AT / £ , 

where (■) denotes averaging with respect to II e . Now, expanding 7r e (including J s ), we derive for any 

<?eLi(n £ ) 

(g) e = (9) - 2 ^((l ^^ du ) (ff))) + • • • • (D-5) 

Hence, using the fact that (<f>) = we obtain that there is a constant C 2 such that 

|E[#y T/e ,z)]|<c a> /E. 

We then conclude from (|D73|) that R £ - 5 is of 0(e 3/2 + eVS + Sy/e). 



E Proof of Proposition 12.51 

Define P*, as the solutions to 



0(1): 


(£ 2 ) Pq.q — 


0, 




*) 


= h(x),' 




(^2> Plfi = 


-V*P * , 




*) 


= 0, 


0(VS) : 


(^2) ^0,1 = 


- (Mi) P * , 


Po,i\T, x, 


*) 


= 0, 


0(e): 


p* 


-^2) 2 P * +P 2 * , 






> 




(£3) ^2*0 = 


-yip * -v*p* , 


Fto(T,x, 


*) 


= 0, 


0(5); 


(£ 2 ) Pq 2 = 


- (Mi) P$ A - M 2 P t , 


^0,2(^1 x > 


*) 


= 0, 


0(V78) : 


(£ 2 > pt, = 


_v*p* 1 -lea 2 p* -(Mi)p* , 


P*,x(T,x, 


*) 


= 0, 



where 

(£*) := <£ 2 ) + VeK 2 D 2 , V* := V - F 2 D 2 . 

These correspond to terms in (|2.4ip of the asymptotic approximation to second order where &(z) is replaced 
by d*(z) defined in (I2.46|) . and the terms containing V 2 are removed. We show that these changes alter the 
accuracy of the approximation only at higher order. 
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First, we note that (Po,o ~ Poo) = ^(v^) smce 

(£ 2 ) (Po,o - P *o) = v^F 2 2) 2 P *o, P ,o(T, a:, z) - P * (T, x, z) = 0. 

Next, we define El' S (t, x, z) by 

E{' S := (Po.o + \/iPi,o + V*Po,i) - ( p o,o + V^Pi.o + V^Poti) , 
the difference in the first order approximations. Note that El' S (T, x, z) = and 

(L 2 ) Ef = [VS(V* + + (Mi)] (P * - Po.o) + eV 2 D 2 Pl + \/^V 2 2) 2 P *i. 

Thus, we conclude that £ , ^ d = 0(e + >/e£). 

Similarly incorporating the order e term, we define E 2 (t,x,y, z) by 

El := (Po.o + y/iPi.0 + eP 2 ,o) - (P *, + V^Plo + ^2*0) • 

From equation (|D.5[) and by using T> 2 (Pq,o — Po 0) = O(v^) one can snow that P|(T, x, ?/, z) = 0(e 3 / 2 ) 
then compute 

(£ 2 ) £| = [(P * + V^Pi.o) - (Po,o + V^Pi.o)] + e-A (Po.o - Po,o) + £ 3/ V 2 D 2 P 2 * . 
Incorporating the order x/eJ term, we define Pf(t, x, z) by 

SI := (P ,i + V?Pi,i) - (PqV + ViPi*i) • 

Note that B|(T, x, 2) = and 

(L 2 ) El = (MO [(P * + VePto) - (Po,o + v^o)] + Vi^ea, (P * - P 0)0 ) + ViV (P * iX - Po.r) 



Now define -Ef (i, x, z) by 



£4 :— Pq,2 — Pq,2- 



Note that Pf (T, x, z) = and 

(£ 2 ) Pf = (Mi) (P *! - P ,i) + M 2 (P * - P , ) + % /iy 2 D 2 P * 2 . 

Finally, 

(£2) (E e 2 + VSE £ 3 + SE{j = (yiV + VS (M^) E{' 5 + e 3/2 V 2 V 2 Pl + ^5V 2 V 2 P* a 

eA + VTs^edA (p * - p , ) + sm 2 (p * - p ,o) 



Hence, we conclude 

E e 2 + Vj£f + <5Pf = 0(e 3/2 + eV^ + y/iS). 
Thus, in (13.31) we can make the following replacements 



a i-> cr* : = +2VJ, V 2 h-> 0. 
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